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ABSTRACT 

Excursion set theory, where density perturbations evolve stochastically with the 
smoothing scale, provides a method for computing the mass function of cosmologi- 
cal structures like dark matter halos, sheets and filaments. The computation of these 
mass functions is mapped into the so-called first-passage time problem in the presence 
of a moving barrier. In this paper we use the path integral formulation of the excursion 
set theory developed recently to analytically solve the first-passage time problem in 
the presence of a generic moving barrier, in particular the barrier corresponding to 
ellipsoidal collapse. We perform the computation for both Gaussian and non-Gaussian 
initial conditions and for a window function which is a top-hat in wavenumber space. 
The expression of the halo mass function for the ellipsoidal collapse barrier and with 
non-Gaussianity is therefore obtained in a fully consistent way and it does not re- 
quire the introduction of any form factor artificially derived from the Press-Schechter 
formalism based on the spherical collapse and usually adopted in the literature. 
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1 INTRODUCTION 



The mass function of dark matter halos is a central object in 
modern cosmology, because of its relevance to the formation 
and evolution of galaxies and clusters. It is therefore impor- 
tant to have accurate theoretical predictions for it, first of all 
when the primordial fluctuations are taken to be Gaussian, 
and then when some level of non-Gaussianity is included. 
Non-Gaussianities are particularly relevant in the high-mass 
end of the power spectrum of perturbations, i.e. on the scale 
of galaxy clusters, since the effect of non-Gaussian (NG) 
fluctuations becomes especially visible on the tail of the 
probability distribution. As a result, both the abundance 
and the clustering properties of very massive halos are sen- 
sitive probes of primordial non-Gaussianities (Matarrese et 
al. 1986; Grinstein & Wise 1986; Lucchin et al. 1988; Moscar- 
dini et al. 1991; Koyama et al. 1999; Matarrese et al. 2000; 
Robinson & Baker 2000; Robinson et al. 2000; Lo Verde et al. 
2008; Maggiore & Riotto 2010c; Lam & Sheth 2009; Gian- 
nantonio & Porciani 2010), and could be detected or signifi- 
cantly constrained by the various planned large-scale galaxy 
surveys, both ground based (such as DES, PanSTARRS and 
LSST) and in space (such as EUCLID and ADEPT) see, e.g. 
Dalai et al. (2008) and Carbone et al. (2008). Furthermore, 
the primordial NG alters the clustering of dark matter ha- 



los inducing a scale-dependent bias on large scales (Dalai 
et al. 2008; Matarrese & Verde 2008; Slosar et al. 2008; Af- 
shordi & Tolley 2008) while even for small primordial NG 
the evolution of perturbations on super-Hubble scales yields 
extra contributions on smaller scales (Bartolo et al. 2005; 
Matarrese & Verde 2009). 

The formation and evolution of dark matter halos is a 
highly complex phenomenon, and a detailed quantitative un- 
derstanding of it can only come through large-scale N-body 
simulations, such as the Millennium simulation (Springel et 
al. 2005). Simulations with non-Gaussian initial conditions 
have also been performed (Grossi et al. 2009; Giannantonio 
& Porciani 2010; Wagner et al. 2010). At the same time, 
some analytic understanding of the process of halo formation 
is also desirable, both for the deeper physical understand- 
ing that analytic models offer, and for their flexibility un- 
der changes of parameters of the cosmological model, shape 
of non-Gaussianities, etc. Analytical derivations of the halo 
mass function are typically based on Press-Schechter (PS) 
theory (Press & Schechter 1974) and its extension (Peacock 
& Heavens 1990; Bond et al. 1991) known as excursion set 
theory (see Zentner (2007) for a recent review). In excursion 
set theory the density perturbation evolves stochastically 
with the smoothing scale, and the problem of computing the 
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probability of halo formation is mapped into the so-called 
first-passage time problem in the presence of a barrier. 

The original formulation of excursion set theory (Bond 
ct al. 1991) makes a number of simplifying assumptions, 
both at the technical level, and concerning the physics of 
halo formation. In particular, at the technical level it is as- 
sumed that the smoothed density field 5 evolves with the 
smoothing scale R (or more precisely with the variance S(R) 
of the smoothed density field) in a Markovian way. How- 
ever, this assumption is correct only if the density field is 
smoothed with a window function which is a top-hat in 
wavenumber space, and with such a smoothing function it 
is difficult to associate a mass M to a region smoothed 
with smoothing parameter R, so in practice it is not pos- 
sible to associate a mass to the dark matter halos identified 
in this way. For any other choice of the window function 
(such as a top-hat in real space, for which the relation be- 
tween the mass M and the smoothing scale R is trivially 
M = (4/3)irR 3 p, where p is the average density of the uni- 
verse), the actual evolution of the smoothed density field 
with R is non-Markovian. At the physical level, the crucial 
simplifying assumption of the original formulation of excur- 
sion set theory is that dark matter halo forms through the 
spherical collapse of initial overdensities. However the ac- 
tual process of halo formation, as revealed by N-body sim- 
ulations, is much more complicated, and involves smooth 
accretion, tidal interactions with the environment, as well 
as violent episodes of collisions with other halos, merging 
and fragmentation. 

In a recent series of papers (Maggiore & Riotto 
2010a,b,c) (hereafter MR1, MR2 and MR3, respectively), 
the original formulation of excursion set theory has been 
extended to deal with the non-Markovian effects which are 
induced either by the use of a realistic filter function, or by 
non-Gaussianities in the primordial density field. The ba- 
sic idea is to reformulate the first-passage time problem in 
the presence of a barrier in terms of the computation of 
a path integral with a boundary (i.e. over a sum over all 
"trajectories" S(S) that always stay below the barrier), and 
then to use standard results from quantum field theory and 
statistical mechanics to express this path integral in terms 
of the connected correlators of the theory. A path-integral 
with boundaries of the kind that we obtain is however not 
a very common object even in quantum field theory or sta- 
tistical mechanics, and in MR1 and MR3 we developed the 
technique for evaluating it perturbatively with respect to 
the non-Markovian and the non-Gaussian effects. This pro- 
vided first of all a rederivation of the results of excursion set 
theory which, from the mathematical point of view, is from 
first principles (for instance the absorbing barrier boundary 
condition, which in the original formulation was imposed by 
hand, comes out automatically in the formalism of MR1). 
Furthermore it allows us to include, at least perturbatively, 
the effect of non-Markovianities and of non-Gaussianities. 
In particular, in MR3 we have shown how to include the ef- 
fect of a non-vanishing bispectrum, while the case of a non- 
vanishing trispectrum was considered in Maggiore & Riotto 
(2010d) (see also D'Amico et al. (2010) for an approach to 
non-Gaussianities which combines our technique with the 
saddle point method developed in Matarrese et al. (2000)). 

Of course this extension of excursion set theory, even if 
it provides an improvement of the original formulation from 



the mathematical point of view, still shares the same phys- 
ical limitations of the original formulations, as long as the 
same model for collapse is used. The model for collapse can 
be improved in different, complementary, ways. A crucial 
step was taken by Sheth, Mo & Tormen (2001) who took 
into account the fact that actual halos are triaxial (Bardeen 
et al. 1986; Bond & Myers 1996) and showed that an ellip- 
soidal collapse model can be implemented, within the excur- 
sion set theory framework, by computing the first-crossing 
rate in the presence of a barrier B C \(S) which depends on S 
("moving barrier"), rather than being constant at the value 
8 C of the spherical collapse, 



Boi(S) 



1 + 0.4 



(!) 



(i) 



Physically this reflects the fact that low-mass halos (which 
corresponds to large S) have larger deviations from spheric- 
ity and significant shear, that opposes collapse. Therefore 
low-mass halos require a higher density to collapse. In con- 
trast, very large halos are more and more spherical, so their 
effective barrier reduces to the one for spherical collapse. 
In order to improve the agreement between the prediction 
from the excursion set theory with an ellipsoidal collapse 
and the N-body simulations, Sheth, Mo & Tormen (2001) 
also found that it was necessary to replace 8 C with ^/a8 c , 
where ^fa ~ 0.84 was obtained by requiring that their mass 
function fits the GIF simulation. The moving barrier there- 
fore becomes 



B cl (S) ~Ja5 c 



1 + 0.4 



\aSV 



(2) 



The parameter a cannot be derived from the dynamics of 
the ellipsoidal collapse. Rather on the contrary, the ellip- 
soidal collapse model would predict a — 1 because in the 
limit S = a 2 — > (i.e. in the large mass limit) halos be- 
come more and more spherical, and therefore the barrier 
must reduce to that of spherical collapse. This mismatch 
might be originated by the fact that, as mentioned above, 
halo collapse is a very complex dynamical phenomenon, and 
modeling it as spherical, or even as ellipsoidal, is a signif- 
icant oversimplification. In addition, the very definition of 
what is a dark matter halo, both in N-body simulations and 
observationally, is a difficult problem. In MR2 it was pro- 
posed that some of the physical complications inherent to 
a realistic description of halo formation can be included in 
the excursion set theory framework, at least at an effective 
level, by taking into account that the critical value for col- 
lapse is itself a stochastic variable, whose scatter reflects a 
number of complicated aspects of the underlying dynam- 
ics (see also Audit et. al. (1997); Lee & Shandarin (1998); 
Sheth, Mo & Tormen (2001) for earlier related ideas). Solv- 
ing the first-passage time problem in the presence of a bar- 
rier which is diffusing around the value S c of the spherical 
collapse model, it was found in MR2 that the exponential 
factor in the Press-Schechter mass function changes from 
cxp{-5 2 /2cr 2 } to exp{-a<5 2 /2cr 2 }, where a = 1/(1 + D B ) 
and Db is the diffusion coefficient of the barrier. The numer- 
ical value of Db , and therefore the corresponding value of a, 
depends among other things on the algorithm used for iden- 
tifying halos. From recent N-body simulations that studied 
the properties of the collapse barrier, a value Db — 0.25 was 
deduced in MR2 predicting a ~ 0.80 (up to a smaller than 
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about 3) We remark that the deduced value of a also holds 
when the collapse is ellipsoidal which was a good fit to the 
average threshold barrier found by N-body data. The value 
of a ~ 0.80 isin excellent agreement with the exponential 
fall off of the mass function found in N-body simulations, 
for the same halo definition. 

The path-integral formulation developed in MR1 and 
MR3 was restricted to the case of a constant barrier 8 C (while 
in MR2 were considered the stochastic fluctuations around 
it). The aim of this paper is to extend the path integral 
formulation of excursion set theory to the case of a generic 
moving barrier, and to provide analytical expressions which 
can be used to calculate the corresponding first-passage time 
probability. 

Given that the Sheth-Tormen (ST) halo mass function 
is widely used in the literature, we believe that it is inter- 
esting to derive it by computing the first-crossing rate with 
an ellipsoidal barrier from first principles. To the best of 
our knowledge, an analytical expression of the first-crossing 
rate was given in Sheth & Tormen (2002) just as a fit to 
the N-body data and its derivation has been sketched only 
recently in Lam & Sheth (2009). As we shall see, this deriva- 
tion is not free from drawbacks. There are other good rea- 
sons why solving analytically for the first-crossing rate with 
a generic moving barrier is interesting. First, excursion set 
theory can be applied to characterize the cosmic web (Shen 
et al. (2006)). Combining models of triaxial collapse with 
excursion set theory, cosmic sheets are defined as objects 
that have collapsed along only one axis, filaments have col- 
lapsed along two axes, and halos are objects in which triaxial 
collapse is complete. Computing the abundances of cosmic 
sheets, filaments and halos within the excursion set theory 
amounts again to solving a first-time passage problem with 
the corresponding moving barriers 



B s hcct(S) ~ \/a S c 

-Bfilam(S') ~ va 5c 



The insertion of each moving barrier into the excursion set 
approach provides estimates of the mass fraction in sheets, 
filaments and halos as a function of mass and time. Secondly, 
moving barriers are adopted in modelling through the excur- 
sion set method the sizes of ionized regions during the epoch 
of reonization (Furlanetto et al. (2004)), while Sheth & Tor- 
men (2002) suggested that moving barriers could effectively 
incapsulate a wide variety of phenomena such as suppres- 
sion of the collapse of small, low-mass, overdense patches in 
models in which dark matter is warm. For a given choice of 
the barrier, the first-crossing rate can in principle be eval- 
uated with numerical techniques (Bond et al. 1991; Zhang 
& Hui 2006), but it interesting to obtain analytic formulas 
valid for a generic functions B(S). Thirdly, as we already 
mentioned, it has become recently clear that detecting a 
significant amount of non-Gaussianity and its shape either 
from the Cosmic Microwave Background (CMB) or from the 
Large Scale Structure (LSS) offers the possibility of open- 
ing a window into the dynamics of the universe during the 
very first stages of its evolution (Bartolo et al. (2004)). It is 
therefore of primary importance to compute the halo mass 



1- 


/ C \ 0.55" 

-°- 56 bf) 
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-°- oi2 (o!) 
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function when NG initial conditions are present. The halo 
mass function with NG has been calculated in Matarrese 
et al. (2000) and Lo Verde et al. (2008) using the PS ap- 
proach with a spherical collapse, while the path integral for- 
mulation of excursion set theory in the presence of NG and 
with a diffusive barrier has been formulated in MR3. The 
main motivation to compute the halo mass function in the 
presence of NG within the excursion set method and with a 
moving ellipsoidal barrier is dictated by the fact that it has 
become customary in the literature to obtain the halo mass 
function with NG by multiplying the ST halo mass function 
with gaussian initial conditions by a form factor obtained 
by dividing the first-crossing rate with NG obtained for the 
PS spherical collapse case (Matarrese et al. (2000) ; Lo Verde 
et al. (2008)) by the PS one (the exception is represented by 
the consistent calculation of MR3, which does not require 
this procedure). It is unclear (at least to us) why and to 
which extent this spurious method should provide a good 
approximation to the correct halo mass function with NG 
and ellipsoidal barrier. The issue is also timely since N-body 
data with NG initial conditions finally exist (Grossi et al. 
2009; Giannantonio & Porciani 2010; Wagner et al. 2010), 
and may be compared to the various theoretical predictions 
for the halo mass functions with NG. They differ at the 
0(20)% level and it is important to understand which error 
is introduced by adopting the form factor procedure. 

The paper is organized as follows. In section 2 we review 
the approach to the computation of the halo mass function 
based on excursion set theory. In particular, in section 2.1 we 
begin with a quick review of the case in which the collapse is 
assumed to be spherical, primordial fluctuations are taken to 
be Gaussian, and the evolution of the density perturbation 
with the smoothing scale is assumed to be Markovian. This 
is the setting considered in the classical paper by Bond et al. 
(1991). We will then proceed toward increasing complexity. 
In Section 2.2 we review the the basic points of the approach 
developed in MR1, MR2 and MR3. In Section 3 we present 
the computation of the first crossing rate for a generic mov- 
ing barrier, while Section 4 contains the generalization of 
the computation to the case of NG initial conditions. Vari- 
ous technical details are collected in Appendices A-D. 



2 THE HALO MASS FUNCTION IN 
EXCURSION SET THEORY 

The halo mass function can be written as 



dn{M) 
dM 



= /(") 



p dlna-^M) 
M2 dlnM 



(5) 



where n(M) is the number density of dark matter halos of 
mass M , a(M) is the variance of the linear density field 
smoothed on a scale R corresponding to a mass M, and p 
is the average density of the universe. The basic problem is 
therefore the computation of the function f(a). 



2.1 Spherical collapse, Gaussian fluctuations, and 
Markovian evolution with the smoothing scale 

Let us summarize the basic points of the original formula- 
tion of excursion set theory. One considers the density field 
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8 smoothed over a radius R, and studies its stochastic evo- 
lution as a function of the smoothing scale R. As it was 
found in the classical paper by Bond et al. (1991), when the 
density S(R) is smoothed with a sharp filter in momentum 
space, and the density fluctuations have Gaussian statistics, 
the smoothed density field satisfies the equation 

dS(S) 



W,S)\ s =s=0. 



(10) 



OS 



= V(S) , 



(6) 



where S = a 2 (R) is the variance of the linear density field 
smoothed on the scale R and computed with a sharp filter 
in momentum space, while r/(S) is a stochastic variable that 
satisfies 



<77(Si)r;(S2)> = MS 1 ! - S 2 ) . 



(7) 



where Sd denotes the Dirac delta function. Equations (6) 
and (7) are the same as a Langevin equation with a Dirac- 
delta noise il(S), with the variance S formally playing the 
role of time. Let us denote by 11(5, S)dS the probability den- 
sity that the variable 8(S) reaches a value between 8 and 
8 + dS by "time" S. A textbook result in statistical physics 
is that, if a variable 8(S) satisfies a Langevin equation with 
a Dirac-delta noise, the probability density 11(5, S) satisfies 
the Fokker-Planck (FP) equation 

^ = i ^H tZ) 

dS 2 dS 2 ' { ' 

The solution of this equation over the whole real axis — oo < 
8 < oo, with the boundary condition that it vanishes at 
8 — ±oo, is 



n°(8,s) 



-5 2 /(2S) 



(9) 



and is nothing but the distribution function of PS theory. 
Since, in hierarchical models of structure formation, as R 
increases, i.e. as the halo mass increases, the variance S de- 
creases monotonically, in Bond et al. (1991) it was realized 
that we are actually interested in the stochastic evolution of 
8 against S only until the "trajectory" crosses for the first 
time the threshold S c for collapse. The threshold value S c is 
estimated within the spherical collapse model where a spher- 
ically symmetric inhomogeneity behaves like a closed col- 
lapsing universe. The underlying idea behind the PS theory 
is that the comoving number density of collapsed haloes can 
computed from the statistical properties of the linear density 
field, assumed to be Gaussian. In this picture haloes form 
when the smoothed linear density contrast is larger than 
S c ~ 1.68 which is obtained computing the linear density 
contrast at the collapse time. This result can be extended 
to arbitrary redshift z by reabsorbing the evolution of the 
variance into 8 C , so that 8 C in the above result is replaced by 
S c (z) — S c (0)/D(z), where D(z) is the linear growth factor. 
Notice that all the subsequent stochastic evolution of 8 as a 
function of S, which in general results in trajectories going 
multiple times above and below the threshold, is irrelevant, 
since it corresponds to smaller-scale structures that will be 
erased and engulfed by the collapse and virialization of the 
halo corresponding to the largest value of R, i.e. the smallest 
value of S, for which the threshold has been crossed. In other 
words, trajectories should be eliminated from further con- 
sideration once they have reached the threshold for the first 
time. In Bond et al. (1991) this is implemented by imposing 
the boundary condition 



The solution of the FP equation with this boundary condi- 
tion is 



n(<5, 5) = 



-« 2 /(2S) 



-{25 c -5) 2 /(2S) 



(11) 



and gives the distribution function of excursion set theory. 
The first term is the PS result, while the second term in 
cq. (11) is an "image" Gaussian centered in 8 — 28 c . Inte- 
grating this 11(8, S) over dS from — oo to S c gives the proba- 
bility that a trajectory, at "time" S, has always been below 
the threshold. Increasing S this integral decreases because 
more and more trajectories cross the threshold for the first 
time, so the probability of first crossing the threshold be- 
tween "time" S and S + dS is given by J r (S)dS, with 

F{S) = ~ f C d8U(S;S). (12) 

With standard manipulations (see e.g. Zentner (2007) or 
MR1) one then finds that the function f(a) which appears 
in eq. (5) is given by 



/(a) = 2a 2 T(a 2 ) 



(13) 



where we wrote S = a 2 . Using eq. (11) one finds the PS 
prediction for the function /(a), 



/ps(o-) = 



5*1/2 



-5 2 /(2<r 2 ) 



,-K/(2S) 



(14) 



Observe that, when computing the first-crossing rate, the 
contribution of the Gaussian centered in 8 = and of the 
image Gaussian in eq. (11) add up, giving the well-known 
factor of two that was missed in the original PS theory. 



2.2 Path integral formulation of excursion set 
theory 

While excursion set theory is quite elegant, and gives a first 
analytic understanding of the halo mass function, it suf- 
fers of two important set of problems. First, it is based on 
the spherical collapse model, which is, as we already men- 
tioned, a significant oversimplification of the actual complex 
dynamics of halo formation. The second set of problems of 
excursion set theory is of a more technical nature, and is 
due to the fact that the Langevin equation with Dirac-delta 
noise, which is at the basis of the whole construction, can 
only be derived if one works with a sharp filter in momentum 
space, and if the fluctuations are Gaussian. However, as it 
is well known (Bond et al. 1991), and as we have discussed 
at length in MR1, with such a filter it is difficult to asso- 
ciate a halo mass to the smoothing scale R. When one uses 
a sharp filter in coordinate space, the evolution of the den- 
sity with the smoothing scale becomes non-Markovian, and 
the corresponding first-passage time problem is technically 
much more difficult. In particular, the distribution function 
11(5, S) no longer satisfies a local differential equation such as 
the FP equation. The issue is particularly relevant when one 
wants to include non-Gaussianities in the formalism, since 
the inclusion of non-Gaussianities renders again the dynam- 
ics non-Markovian. Neglecting the non-Markovian dynamics 
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due to the filter function would lead to incorrectly assigning 
to non-Gaussianities in the primordial density field effects 
which are rather due, more trivially, to the procedure that 
one has adopted for smoothing the density field. 

In MR1,MR3 has been developed a formalism that al- 
lows us to generalize excursion set theory to the case of 
a non-Markovian dynamics, either generated by the filter 
function or by primordial non-Gaussianities. The basic idea 
is the following. Rather than trying to derive a simple, local, 
differential equation for 11(5, S) (which, as shown in MR1, 
is impossible; in the non-Markovian case 11(5, S) rather sat- 
isfies a very complicated equation which is non-local with 
respect to "time" S), we construct the probability distribu- 
tion 11(5, S) directly by summing over all paths that never 
exceeded the threshold 5 C , i.e. by writing 11(5, S) as a path 
integral with boundaries. To obtain such a representation, 
we consider an ensemble of trajectories all starting at So = 
from an initial position 5(0) = 5o and we follow them for a 
"time" S. We discretize the interval [0, S] in steps AS — e, 
so Sk ~ ke with k = 1, . . . n, and S„ = S. A trajectory is 
then defined by the collection of values {S 1 , . . . , 5 n }, such 
that 5(Sk) = 5k- The probability density in the space of 
trajectories is 



W(5 ;5i 



, S„; S n ) = (Sd(6(Si)-6i) . . . S D (5(S n )-S„)) ,(15) 



where Sd denotes the Dirac delta. Then the probability of 
arriving in 5 n in a "time" S„, starting from an initial value 
5o, without ever going above the threshold, is 

n e (5o;5 n ;S'„) = / dSi . . . / dS n -i 



xW / (5 ;5i 



, 5 n -l, 8 n ; S n )- 



(16) 



The label e in IT reminds us that this quantity is de- 
fined with a finite spacing e, and we are finally inter- 
ested in the continuum limit e — > 0. As discussed in 
MR1 and MR3 (see Eqs. (23)- (27) and discussion therein), 
W(8 ; 5i, . . . , 5„_i, 5„; S n ) can be expressed in terms of the 
connected correlators of the theory, 



(17) 



W(6o;6i,...,6 n ;S n )= \ VXe 

where 
V\\ 

and 

n 

Z — i \^ Xi8 



d\i d\„ 

~2Y '" "2V 



(18) 



(19) 



We also used the notation Si — 5(Si), and (5i . . . S n )c denotes 
the connected n-point correlator. So 



II e (5o;5„; S n ) 



d5i . . . d5n-i / T>\ 



(20) 



1 In eqs. (9) and (11) wc had implicitly assumed 8q = 0. In the 
following however it will be necessary to keep track also of the 
initial position <5q. 



When 6(S) satisfies eqs. (6) and (7) (which is the case for 
sharp filter in wavenumber space) the two-point function 
can be easily computed, and is given by 



(5(S i )S(S j ))=imn(Si,S j 



(21) 



If furthermore we consider Gaussian fluctuations, all n-point 
connected correlators with n ^ 3 vanish, and the probability 
density W can be computed explicitly, 



W gm (<5o;5!, 



■s)- 1 e -*E:->^ 

,i,n) ~ (27T6)"/2 



Si)- 



(22) 



where the superscript "gm" (Gaussian-Markovian) reminds 
us that this value of W is computed for Gaussian fluctua- 
tions, and when the evolution with respect to the smooth- 
ing scale is Markovian. Using this result, in MR1 we have 
shown that, in the continuum limit, the distribution func- 
tion II e =o(<5; S), computed with a sharp filter in wavenumber 
space, satisfies a Fokker-Planck equation with the boundary 
condition n e =o(<5 c , S) = 0, and we have therefore recovered, 
from a path integral approach, the distribution function of 
excursion set theory, eq. (11). Considering a more realis- 
tic filter, such as a step function in coordinate space, nec- 
essarily introduces non-Markovianity and the computation, 
which is quite non-trivial from a technical point of view, has 
been discussed in great detail in MR1. In order to make the 
computation of the first-crossing rate with a moving barrier 
more clear, from now on we will adopt the step function 
in wavenumber space as a filter and eliminate the source of 
non-Markovianity given by the choice of the window func- 
tion. The effect of a more realistic filter function could then 
be computed as in MR1. The effect, however, will be tiny 
and totally negligible in the large mass range we are mostly 
interested in for the non-Gaussian case. Let us just close 
this subsection by reminding the reader about some useful 
properties of the path integral formulation which will turn 
out to be useful in the following. We will encounter objects 
such as 

Y,F(Si) dS 1 ...dS„-id l W gln (S ;S 1 ,...,Sn;Sn), (23) 

«=i ^-°° 

where F denotes a generic function. To compute this expres- 
sion we integrate di by parts, 

dSl... dSn-l cW gm (5 ; 5l, . . . , Sn, Sn) 

' —00 

r s c 

dSi...dSi... dSn-i (24) 

xW(5o;5i,...,Si = Sc,- ..,S n -i,S„;S n ) , 

where the notation dSi means that we must omit dSi from 
the list of integration variables. We next observe that W sin 
satisfies 



iy gm (5 ;5 1 ,...,5 1 = 5 c ,...,5„;5„) 

H/ gm (5 ;5 1 ,...,5,_ 1 ,5 c ;^) 

xW sm {S c ;S i+1 ,...,Sn;Sn-S l ), 



(25) 



as can be verified directly from its explicit expression (22). 
Then 

Sc fSc 

dSi...dSi-i / d5i+i . . . d5n-i 
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xW sm (S ; Si,..., Si-i, 5 C ; S t )W gm (S c ; S l+1 , ...,6 n ;S n - Si) 
= TJf n (So;S c ;S i )Uf n (Sc;S n ;S n - Si) , (26) 

and to compute the expression given in eq. (23) we must 
compute objects such as 



n-l 

E 



F(Si)Ur (So; S c ; &)nf m (<5 c ; S n ; S n - Si). 



(27) 



To proceed further, we need to know IIf m (<5o; S c ; Si). By def- 
inition, for e = this quantity vanishes, since its second 
argument is equal to the the threshold value 5 C , compare 
with eq. (10). However, in the continuum limit the sum over 
i becomes 1/e times an integral over an intermediate time 
variable Si, 



n_1 i rS n 



dS t , 



(28) 



so we need to know how ITf m (<5o; S c ; Si) approaches zero 
when e — > 0. In MR1 we proved that it vanishes as \fi, 
and that 

nf m (5 ; S c ; S) = ^^= S -^l e -(^of/(2s ) + 0{e) (29) 

Similarly, for 5„ < 8 C , 

nr(8 c ; <5„; S) = v^ 4= Ss^£L e -«c-s»)>/(?s) + Q{e) (30) 

In the following, we will also need the expression for IIf m 
with the first and second argument both equal to S c , which 
is given by (see again MR1) 



nr(8c;8 c ;S) = 



2-kS 3 ' 2 



(31) 



The two factors y/e from cqs. (29) and (30) produce just an 
overall factor of e that compensates the factor 1/e in eq. (28), 
and we are left with a finite integral over dSi . Terms with two 
or more derivative, e.g. didj, or di,djdk acting on W, with all 
indices i,j, k mailer than n, can be computed similarly, and 
have been discussed in detail in MR1. With these technical 
details in mind, one can proceed to the computation of the 
first-crossing rate in the presence of a moving barrier. 



3 PATH INTEGRAL WITH MOVING 

BARRIER: GAUSSIAN FLUCTUATIONS 
AND MARKOVIAN EVOLUTION WITH 
THE SMOOTHING SCALE 

In this section we discuss the first-crossing rate for a generic 
moving barrier B(S), specializing to the ellipsoidal one at 
the end. We consider first the case of Gaussian primordial 
fluctuations, and we will assume that the evolution with 
the smoothing scale is Markovian. Similarly to the constant 
barrier case, the probability of arriving at 8 n in a "time" S„, 
starting from the initial value So — 0, without ever going 
above the threshold, is 



IT e (<5„; S n 



B(Si) /-B(S n _i) 

dSi . . . / d8 n -i (32) 

■oo J — oo 

xW(5 ;Si,. . . ,S„-i,S„;S„). 



Since we are considering the Gaussian and Markovian case, 
W(8o] Si, ... , S„-i, S n ; S n ) can be expressed in terms of the 
connected two-point function of the theory, as 



VX 



W{5 ;Si,...,5 n ;S n )-. 

x exp < i \J XiSi — - \J XiXj min(S'i, Sj) 



(33) 



ij=i 



Taking the derivative with respect to the time S n = S of 
eq. (32) and using the fact that iXj (j — 1, • • • , n) can be 
replaced dj, we discover that n e (<5;S) satisfies the Fokker- 
Planck (FP) equation 



dU e (5;S) _ ld 2 U e (5;S) 
' 2 



i)S 



<)<)- 



(34) 



(where we used the notation S n — S). To determine the 
boundary condition to be imposed on the solution of eq. (34) 
we proceed as follows. We start from eq. (32), with W given 



by eq. (22) and, shifting the variables Si (i 

Si — > Si — B(Si), we obtain 



, n) as 



n € {8n + Bn\ Sn) = / dSi . 

J — oc 

i 



dSn-l 



(27re)™/ 2 

(-0 

dSn-l- 

' — oo 



'2-Kt 

xll e (5 n _i + B n -i;S n -i) 



e -^-[S n -S n _ 1+ B„ 



(35) 

Now 



where we used the notation Bi = B(Si), so B„ = B n 
let S„_i = S so Sn - S + e, and 5 n + B(S) = S,S n - 
8n-i — AS. For fixed 5 n , we have dS n -i = — d(A5). By 
further taking the limit e — > (assuming that -6(5*) is a 
continuous and difierentiable function), eq. (35) becomes 

foe 

n € =o(8;S)= / d(AS)S D (AS)U e=0 (S-AS;S). (36) 

JS-B(S) 

From this relation we get the boundary condition. If S — 
B(S) the integral is over half of the support of the Dirac 
delta and so H e=0 (B(S); S) = (l/2)U e=0 (B(S);S) hence 
U e=0 (B(S); S) = 0. Furthermore, if S > B(S), the support 
of the Dirac delta is outside the integration limits and there- 
fore we conclude that 



n, =0 (S;S) = for S^B(S). 



(37) 



In the continuum limit the first-crossing rate is then given 
by 



T(S) = 



d_ 

' dS 

dB(S) 
dS 



HS) 



d5U^o{S;S) 



U e=0 (B(S),S) 



(38) 



B(S) 



<IS 



dn t = (8;S) 
ds 



The first term on the right-hand side vanishes because of the 
boundary condition, while the second term can be written 
in a more convenient form using the FP equation (34), so 



ns) 



B(S) 



,18 



d 2 U e=0 (S;S) 
OS 2 
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l dU e=0 (6;S) 

2 



05 



(39) 



6=B(S) 



To compute the probability n e= o(<5 n , S n ) we proceed in the 
following way. At every i-th step of the path integral we 
Taylor expand the barrier around its final value 



B 



(p) 



B(Si) = B(S n ) + J2 -T~ (Si- S n ) 



p=i 



where 
B 



(p) _ d?B(S n ) 
' dS p 



(40) 



(41) 



(so in particular B 



«>j 



B(S„ 



shift in the variable Si (i — 1, 



= Bn). We now perform a 
, n — 1) in the path integral 



<5; — > <5; 



^ p! 

P =i 



^Oi OnJ 



Then II e (5 n ; S n ) can be written as 



U e (5„; S n 
where 
Z = i 



d5i. 



dSn-i / T)\ e 



(42) 



(43) 



1 = 1 



n 



».j'=i 



n — 1 oc (p) 

+*E A »E^^-^) P - 



(44) 



i — 1 p — 1 

We next expand 




\Ji On) 






i+iE^E^-^-*')* 



(45) 



i — 1 p — 1 



?(p) d(<?) 



p,q=l 



('Ji ^nj (^j ^nj i ' 



"2L XiXj 

and we write n e (5„; SV,) as 

n e (5 n ;S„) = n™( ( 5 n ;S*„)+n( 1) ( ( 5 n ;S'„) 

+n( 2) (5 n ;S' n ) + --- , 

where 



(46) 



r(0) 



U^ (S n ;S n ) = 



He (^ni »-5ti 



and 



\/2nS n 



-**/(2S n ) 



-(2S„-<5„) 2 /(2S„) 



We have therefore formally expanded H e= o(8„,S„) in a se- 
ries of terms IT^ =0 , II^_ , etc., in which each term is itself 
given by an infinite sum over indices p, q, .... To proceed fur- 
ther, we must either perform some approximation, or iden- 
tify a suitable small parameter, and organize the terms in 
a systematic expansion in such a small parameter. In the 
next subsections we first discuss the approximation in which 
one can rederive the Sheth-Tormen result, and we will then 
compare it with two complementary, and more systematic, 
expansions. 



3.1 The Sheth-Tormen approximation 

To attack the problem, a first idea is to perform the integrals 
in eqs. (48) and (49) approximating (S n — S , i ) p_1 ~ S^ 1 
inside the integrals. This is in fact equivalent to the approx- 
imation made in Lam & Sheth (2009), see in particular their 
eq. (20). The detailed calculations, within our formalism, are 
reported in Appendix A and one obtains the first-crossing 
rate for a moving barrier 



.Fst(S) 



,-B 2 (S)/(2S) 



27rS 3/2 Z^ p \ 
p=0 



(_Sy d p B(S) 



OSp 



(50) 



This expression agrees with the one suggested in Sheth & 
Tormcn (2002). Notice that for the cases of constant barrier 
B(S) = 8 C and of a linear barrier B(S) = 8 C + /3S, which 
are the known examples where the first-crossing rate can be 
computed analytically by solving exactly the FP equation in 
the presence of such a barrier (for the linear barrier see Sheth 
(1998) and Section IX of Zentner (2007)) the first-crossing 
rate (50) reproduces the correct answer. When applied to 
the ellipsoidal barrier given in eq. (2), and restricting the 
sum to p ^ 5, one recovers the ellipsoidal collapse result of 
Sheth & Tormen (2002) 



Ht(S) 



VaSc -B 2 (S)/(2S) 



2ttS 3 / 2 
5 



1 + 



-B-OQ' 



fa 5c 



p=0 

_ J -B 2 (S)/(2S) 



2nS 3 / 2 



0.067 



aS'i 



(51) 



B 



(p) 



This procedure is, however, not free from drawbacks. In- 
deed, the restriction of the sum to p ^ 5 is not justified and 
is merely dictated by the fact that stopping arbitrarily the 
series at p = 5 provides a a good fit to the N-body simula- 
tions. 2 However, if the sum over p is extended up to infinity 
(47) the sum simply resums to -8(0) since, performing a Taylor 
expansion of B(So — S) in powers of S and setting finally 
S*o = S, we have 



E/ d *-"<«»-i2J!r (48) 

«=i J -° c P =i 

x (Si - S n ) p cW gm (<5 ; 61,..., S n ; S n ) . 



5(0) = 



p=0 



-S) p d p B(S) 
p\ dSP 



Since B(0) = \/aS c , we just end up with 



u'r' (.*,>; s n ) = I J2 



tWi 



dSi . . . d5„ 



E 



B (p) B (i) 



•S p— OO \& ) 



\Ja8 c 



-B'(S)/(2S) 



i,j=l" °" p,q=l 

x (S t - S n ) p (S, - S n ) q d i d j W em (8 ; Si, . . . , 6 n ; S n ) . (49) 



2 We thank Ravi Sheth for discussions about this point. 



(52) 



(53) 
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so the correction term ~ S in eq. (51) seems an artifact of 
stopping the sum to p — 5. This is a rather puzzling result, 
since this correction is known to fit well the data, and is 
widely used in the literature. This calls for a different and 
more rigorous approach where the integrals are performed 
without the approximation (S n — S'i) p_1 ~ S^ -1 . We discuss 
two different possible approaches in the next two subsection. 



3.2 Expansion of II e (<5, S) in derivatives of B(S) 

In order to develop a more systematic expansion, we first 
consider the case of a barrier B(S) which is slowly varying 
with S. In this case, the small parameters are the derivatives 
of the function B(S). 

At first one might think that such an approximation, 
altough useful in some cases, would not apply to the barrier 
which corresponds to the the ellipsoidal collapse, eq. (2). In 
this case infact B e \(S) is given by a constant plus a term pro- 
portional to S J with 7 ~ 0.6 < 1, and therefore already its 
first derivative, which is proportional to S* 7-1 is large at suf- 
ficiently small S, and formally even diverges as S — > 0. How- 
ever one should not forget that, in practice, even the largest 
galaxy clusters than one finds in observations, as well as in 
large-scale JV-body simulations, have typical masses smaller 
than about 10 15 /i _1 M Q which, in the standard ACDM cos- 
mology, corresponds to values of S = <t 2 (M)>0.35, see e.g. 
Fig. 1 of Zentner (2007). Even for such a value, which is the 
smallest in which we are interested, the value of B' cl (S) is 
just of order 0.3 which means that, in the range of masses of 
interest, the barrier of ellipsoidal collapse can be considered 
as slowly varying. 

We therefore expand H e (5 n ; S n ) in powers of the deriva- 
tives of the barrier, keeping terms with the same number of 
derivatives, so for instance a term proportional to d 2 B/dS 2 
is taken to be of the same order as (dB/dS) 2 . Working up 
to terms of second order in the derivatives we get 

n e (s n ;S n ) = ni 0) (s n ;S n ) + n ( c a) (8 n ;S n ) + ni b) (s n -s n ) 



+n[ c) (s n; s n ), 

where 

n-l 

ni a \s n ;S n ) = ^KOSi-Sn) 

x / dSi . . . d<5 n _ 1( W gm . 

J — oo 
n-l 

nf\8 n ;S n ) = i J^B^ -5„) 2 

i=l 

x / dSi . . . dS n -idiW sin , 



(54) 



(55) 



(56) 



n-l 



n<f\5 n ;S n ) = ±J2 (K)\Si-S n )(Si-S n ) (57) 



r-B„ 



and we used a prime to denote the derivatives of B(S n ) with 
respect to S„. Observe that IT"' and iT 6 - 1 are linear in the 
first and second derivative, respectively, and come from the 



terms p = 1,2 of IT 1 ', while IT C ' is quadratic in the first 
derivative, and is the term p = q = 1 of IT 2 ' . 

In Appendix B we compute these three terms, in the 
continuum limit, using the techniques developed in MR1. 
For the first term we find 



n?J (S n ;S n ) 



-25; 



/ (B n — S n ) 

\/2nS n 



-(2S„-«„) 2 /(2S„) 



(58) 



Observe that it satisfies the boundary condition 
n^ ((5 n ;5 n ) = when 8 n = B n , as it should. For 
the second term we get 



n^ (5 n ;S n ) 



v2ttS„ 



1 



B n (B n — S n ) 



2tt 

(2S„-i„) 2 /(2S„) 



(59) 



7r_B n Erfc 



/ 2B n — S ri 



and again vanishes linearly as S n — > B n - The third term is 
given by 

,/\2 (Bn-Sn) r -(2B n -S n ) 2 /(2S n ) l 



,; ,(60) 



ni%(5n;Sn) = ~2(B'ny 

\/2ni>n 

and vanishes quadratically as 8 n — > B n - This means that 
in the end it does not contribute to the first-crossing rate, 
since, using eq. (39), the latter is given by the derivative of 
II e= o(<5 n ; S n ) with respect to S n , evaluated in S n — B n . 

It is interesting to check explicitly that this solution for 
n(8 n ;Sn) satisfies the FP equation, up to order to which 
we have computed, i.e. up to terms of second order in the 
derivatives of the barrier, included. Define the FP operator 

l#_ 

2 381 ' 



D 



d 

8Sn 



(61) 



and define / (0) ,.../ (c) from 



DUf =0 (8 n ;S n ) 



TV q 



3/2' 



i -(2B n -S) 2 /(2S n )rA 



(62) 



where A = (0), (a), (b), (c) so, up to terms of second order 
(included) in the derivatives of the barrier, 



DU e=0 (8 n ;S„) 



* S: 

f(0) 



3/2 



= -(2B„-5„) 2 /(2S„) 



x[.r> + r> + .r + f 



(b) 



rWi 



(63) 



Inserting the expressions for II 
puted above we get 



rrw Yt<,"> rr (c) com- 



(0) 



r(») 



,('') 



/( 0) 
/( c) 



= (2B n — 8„)B n , 

= —(2B n — 8 n )B n — S n (B n — 8 n )B n 

+ [2(B n - 8n)(2B n - Sn) - S n ](B n ) 2 
= Sn(B n — 8n)B„ + 0(B n , B n B n , (B. 
= -[2{B n - Sn)(2B n - Sn) - S n ](B, 

+£ > (B n ,B n B n ,(B n ) ) 



n) 3 ) 
/ n2 



(64) 

(65) 
(66) 

(67) 

Therefore the sum IT^ + II^,, + II^ + Il[% satisfies the 
FP equation, modulo terms of third order in the derivative 
of the barrier. 

The first-crossing rate is then readily evaluated through 
eq. (39). The zero-th order contribution from II^_ is 



T (U >(S) = 



BjS) 

2nS 3 / 2 ' 



-B 2 (S)/(2S) 



(68) 
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Figure 1. The Sheth-Tormen first-crossing rate for the ellipsoidal 
barrier J^S^ (dashed black line), compared to the first-crossing 

(2) 
rate J-j (solid blue line) obtained from the expansion in deriva- 
tives of the barrier, as a function of v. 



while the higher orders give 



T W {S) 



-B / ('5') c -a(S) 2 /(2S) 



4tt 



-B(S) 2 /(2S) 



7rB(5)Erfc 



and J r ' c ' = 0, as already mentioned. In Fig. 1 we compare 
the Sheth-Tormen first crossing rate J-st(S) to the quantity 

(2) (5) = T m (S) + T (a) (S) + T (b) (S) , (71) 



?: 



i.e. to the first crossing rate obtained by performing the 
expansion in derivatives of the barrier, up to second order 
(included) in the derivatives, while in Fig. 2 we plot the rela- 
tive difference (J^J- -Fst)/-^*St. We see that the two results 
agree perfectly at large values of v (i.e. at large masses), and 
they still agree to better than 10% down to v = 1. 

(2) 

The fact that the J-^ cr is numerically quite close to J-st 
provides a more satisfying derivation of the ST mass func- 
tion, showing that the approximation (S„ — 5'i) p_1 ~ S^ -1 , 
together with the truncation top = 5 of the series in eq. (50), 
in the end gives a simple analytic formula which is numer- 
ically quite close to the result of a derivation based on a 
systematic expansion. 

For comparison, we also report in Figure 3 the first- 
crossing rate for filaments (blue), sheets (red) and halos 
(brown). The dashed lines refer to the ST approximation 
(50) with p ^ 5, while the continuous ones refer to our re- 
sult (71). 



3.3 Expansion of II e =o(5n, S n ) in powers of (B n — S„) 

In this subsection we describe a different expansion scheme, 
which allows us to resum a large number of terms. The basic 
idea is that, even if the computation of the distribution func- 
tion II can be interesting by itself in a more general context 
(since the probability distribution of a random walk in the 
presence of a moving barrier is a problem interesting in its 
own right in statistical physics) , for the computation of the 
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Figure 2. The ratio (.Fl — -FstV-^ST, as a function of v. 
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Figure 3. First-crossing rate for filaments (blue), sheets (red) 
and halos (brown). The dotted lines refer to the ST approximation 
(50) with p Sj 5, while the continuous ones refer to our result (71). 



halo mass function we are really interested only in the first- 
crossing rate. Then eq. (39) shows that, in the Gaussian and 
Markovian case, we only need the derivative dTl/dS n evalu- 
ated at 5 n — B n - As shown in eq. (37), Ii(S n ,S n ) vanishes 
at 5 n — B n , so its Taylor expansion around S n — B n starts 
from a term linear in (5 n — B), followed by terms of order 
(5 n — B) 2 , etc. When we compute dW/d8 n in 8 n = B n , the 
terms quadratic and higher-order in (5 n — B), give zero, so 
we do not need the full function II, but only the term linear 
in (5 n — B) in its Taylor expansion around S n — B„. This 
simplifies our task considerably. 

We first compute the part linear in (S n — B n ) of TiS ' . 
Using the results of the previous section, in particular eqs. 
(24), (29) and (30), nP(d n ;S n ) can be rewritten as 



r(l) 



n£ (*„;S , „) = 



B n (B n -s n) y(-rf B() 

p=l 
Jo 



s 



3/2 



xe -BZ/(2S i ) e -(.B n .-S n y/l2(S n -S i )] 



For p — 0, 1 this integral can be computed analytically, see 
appendix C, but for p Js 2 we have not been able to com- 
pute it exactly. However, for our purposes it is sufficient to 
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(i) 



observe that in this expression for n^ =0 there is already a 
factor (B n — 5 n ) in front of the integral over dSi, and the 
integral converges at Si — S n for all p ^ 1, even if in the 
integrand we set 5 n = B n . Therefore 



rd) 



n^ (<5„;S„) 



B n {B n - S n ) v- (-l) p 



E 
P =i 



(?) 



dSi 



(Sn 



B 



\p-(3/2) 



(73) 



-B*/(2S<) 



o S, 

2 



3/2 



+o(B n - s n y 

In appendix C we show that for p = 1 this integral is ele- 
mentary while for p ^ 2 it can be computed in terms of the 
confluent hypergeometric function U(a,b,z). As a result, 



r(l) 



n^ (<5„;S„) 



V2 B„ — 5„ 



,1/2 



-B*/(as„) 



E 

Lp=i 



.1)p 



B^sr'rip-^jyjp-i 



i i£ 

2'2S n 



+C(B„ - 5 n f 



(74) 



where the term p = 1 can be written in a more elementary 
form using £/(0, 6, 2) = 1 and F(l/2) = \pK. Along the same 
lines, we have also computed the generic m-th order (m ^ 1) 
of the expansion of n e =o (see App. D), at the linear order 
in B n — 8 n , and it is given by 



n 



(mi 



(B n -5n)e 35 

«™_1 3 ~" 

m! 2 2 'i 2 



E 

Plv,Pm=l 



B 



(pi) 



■B 



(p m ) 



,S. 



(_i)Er=i pfc+m+i 



E pfc ~ 




m + l 1 B n 

Zi _ Znjf-, 



+ 0{B„- Sn 

where the coefficients c v 



(75) 



l Pj9i ... can be computed by the re- 
cursion relations (DlO)-(Dll). This expression is useful for 
numerical evaluation, but not very illuminating from an an- 
alytic point of view. So it can be useful to keep in mind 
that in the limit 25,! <C B n , i.e. for large halo masses, the 
confluent hypergeometric U function simplifies to 



U k 



1 K 



(2Sn\ 

V Bl ) 



" 1 + °(t 



(76) 



n 



(m) 



We 



The total probability is given by II = ^ c _ n 
have not been able to resum all the terms of the expan- 
sion, but the first few terms are sufficient for the first- 
crossing rate. In fact, the first-crossing rate is readily evalu- 
ated through eq. (39). The zero-th order contribution from 
Ilj, =0 is given by eq. (68) while higher-order contributions 
T m ' are obtained from II^™ () in eq. (75), and are easily 
evaluated numerically. In Fig. 4, we plot _F' ' (blue) and 
J-^°' +J- 1 ' +J- 1 ' 1 + ■ • ■ (red), for the ellipsoidal barrier given 
in eq. (2). We deduce that the sum for II converges quickly 
and the terms after the second one contribute negligibly to 
the first-crossing rate. It is therefore an excellent approxima- 
tion to consider the first-crossing rate for a generic moving 




Figure 4. First-crossing rate for the ellipsoidal barrier (2). ^st 

(dashed black), J^f ' (solid blue), J 7 ^' + J^ 1 ' + J^ 2 ' H (solid 

red). The spherical collapse model, with the same value of a = 
0.707, corresponds to the dotted black line. 



barrier B(S) as given by J r ' ' + -F' , i.e. 

-S 2 (S)/(2S) 
v/2^5' 3 / 2 



HS) = 



B(S) 



+ 



(-S) p d p B(S) r (p - 



p=l 



P! 



OSp 



V 



1 K 



(77) 



For comparison, we also report in Fig. 4 the first-crossing 
rate of the spherical collapse model (dotted line) and the 
Sheth & Tormen (2002) result of eq. (51) (dashed line). Note 
also that eq. (77) reproduces the exact known results for 
the cases of constant and linear barrier shapes. It is also 
interesting to note that J-(S) in eq. (77) and the rate J-^ CI (S) 
computed in the previous section differ by less than 5% for 
v ^ 0.2, for the ellipsoidal barrier (2). It is then reassuring to 
see that our two approaches to the computation of the first- 
crossing rate lead to consistent results, and their difference 
allows us to get a quantitative idea of the theoretical error in 
the computation. The fact that both results are numerically 
quite close to the ST mass function also provides a more 
satisfying justification of the ST mass function itself. 

Armed with these results, we may now proceed to eval- 
uate the halo mass function in the case in which non- 
Gaussianity (NG) is present. 



4 THE ELLIPSOIDAL COLLAPSE AND 
NON-GAUSSIANITY 

Deviations from Gaussianity are encoded, e.g., in the con- 
nected three- and four-point correlation functions which are 
dubbed the bispectrum and the trispectrum, respectively. 
A phenomenological way of parametrizing the level of NG 
is to expand the fully non-linear primordial Bardeen grav- 
itational potential $ in powers of the linear gravitational 
potential <I>l 



$ = $L + /ML ($L - (*L>) 



(78) 



The dimensionless quantity /nl sets the magnitude of the 
three-point correlation function (Bartolo et al. (2004)). If 
the process generating the primordial NG is local in space, 
the parameter /nl in Fourier space is independent of the 
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momenta entering the corresponding correlation functions; 
if instead the process which generates the primordial cos- 
mological perturbations is non-local in space, like in models 
of inflation with non-canonical kinetic terms, /nl acquires 
a dependence on the momenta. The strongest current limits 
on the strength of local NG set the /nl parameter to be in 
the range —4 < /nl < 80 at 95% confidence level (Smith et 
al 2010). 

In MR3 the effect of primordial NG on the halo mass 
function was computed, using excursion set theory, for the 
case of a spherical collapse with constant barrier. In the pres- 
ence of NG the stochastic evolution of the smoothed density 
field, as a function of the smoothing scale, is non-Markovian 
and beside "local" terms that generalize Press-Schechter 
(PS) theory, there are also "memory" terms, whose effect on 
the mass function have been computed using the formalism 
developed in MR1. When computing the effect of the three- 
point correlator on the mass function, a PS-like approach 
which consists in neglecting the cloud-in-cloud problem and 
in multiplying the final result by a fudge factor ~ 2, is in 
principle not justified. Indeed, when computed correctly in 
the framework of excursion set theory, the "local" contribu- 
tion vanishes (for all odd-point correlators the contribution 
of the image Gaussian cancels the Press-Schechtcr contribu- 
tion rather than adding up), and the result comes entirely 
from non-trivial memory terms which are absent in PS the- 
ory. However it turns out that, in the limit of large halo 
masses, where the effect of non-Gaussianity is more rele- 
vant, these memory terms give a contribution which is the 
the same as that computed naively with PS theory, plus sub- 
leading terms depending on derivatives of the three-point 
correlator. 

The goal of this section is to compute, using excursion 
set theory, the halo mass function in the presence of NG and 
for the ellipsoidal collapse, thus extending the findings of 
MR3 obtained for the spherical collapse. This computation 
is motivated by the fact that in the literature the halo mass 
function for the more realistic case of the ellipsoidal collapse 
is obtained, when NG is present, by multiplying the first- 
crossing rate (51) by a form factor 72.(/nl,S') obtained by 
dividing the first-crossing rates with and without NG for the 
PS spherical collapse case 

?ST(fsL,S) = ^st(/nl = 0,5')^(/nl,5') 

= ■Fst(/nl = 0,A) — — . (79) 

>ps(/nl = u,b) 

This procedure has however no rigourous justification and 
its validity should be tested with an explicit computation. 

Similarly to the Gaussian case, the probability of ar- 
riving in S„ in a "time" S n , starting from the initial value 
So = 0, without ever going above the threshold, in the pres- 
ence of NG is given by 



B(Si) 



B(S„-l) 



U. e (d n ;S„) = I dSi . . . / dS„-i 

J — oo J — oo 

xW NG {So;di,. . . ,S n -i,6 n ;S n ). (80) 

where 



W / ng(5o;5i 



j Vn'i &n) 



vx 



I n n 

exp I i 2_^ KSi — - 2_^ A»Aj min(£i, Sj 



»ij'=i 



x exp ■ 



H) 3 
6 



/ {8iSjSk} c Xi\j\k 



(81) 



i,j,fc = l 



We now perform the shift (42) in the Si (i — 1, • • • , n — 1) 
variables and expand the NG contribution to first order 

IT (S n ; s„ ) = n^ {Sn ; s n ) + ni% (S n ;S n ) 
+ n^2 (8n; S n ) + ■ ■ ■ 



dSi ... / dS„-i 

— oo J —oo 



i,j,fc = l 



(SiSjS^cdidjdkWmbiSo; (Si, ... , 5„-i,6 n ; S n ), 

(82) 

where Wmb is the probability density in the space of trajec- 
tories with a moving barrier, so that 



dSi 



B„ 



d8n-i W m b(5o;Si 



, O n \ On) 



TT(0) , tt(1) , n (2) , 



(83) 



In principle the contribution from NG can be computed sep- 
arating the various contributions to the sum according to 
whether an index is equal or smaller than n. In this way, 
however, the computations faces some technical difficulties. 
Fortunately, as discussed in MR3, the problem simplifies 
considerably in the limit of large halo masses, which is just 
the physically interesting limit. Large masses mean small 
values of S„. The arguments Si, Sj and Sk in the correlator 
(SiSjSk) = (S(Si)S(Sj)8(Sk)} c range over the interval [0, 5 n ] 
and, if S„ goes to zero, we can expand the correlator in a 
multiple Taylor series around the point Si = Sj = Sk — S n . 
We introduce the notation 



G 



(P,1,r) 



(S n 





L i 3 fc 


Si = Sj =. 


Then 




W 


S*,)5(5 3 )5(5 fc )) -- 


oo 

= E ^ 

p,q,r= 





(84) 



(S n - S % ) p (85) 



X {Sn Sj) {S n — Sk) CjT3 ' ' (*->ti 



The leading contribution to the halo mass function is given 
by the term in eq. (85) with p — q — r — and we neglect 
subleading contributions, which can be computed with the 
same technique developed in MR3. The discrete sum reduces 
to (S n ) c X^™ • fe-i didjdk and we can split it as 



22 didjd k — d n + 3 



i,j,fc = l 



n— 1 n— 1 

} j dtdjdn + 3 Yj did„ 

j,j — 1 i=l 



n-1 

+ 2_^ didjdk . 

i,j,fc = l 



(86) 
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When applying these derivatives to the W m b, one can use 
the identities proven in MR1 and MR3, namely 



E 



dSi . . . d5 n -i diWmb 



d 
dB, 



il =0 . 



" — 1 [B n 

2J / ddi . . . dd„-i didjWmb — 

iij=1 J-oo 

and 

"-1 rB n 

} j / ddi . ..dSn-ididjdkWmb 



d 2 

dB 2 (S n 



■n e=0 . 



(87) 



(88) 



<f> 



dB 3 (S n 



n £=0 . (89) 



The probability density (82) calculated in this way vanishes 
at the barrier point 8 n = B„ , when one properly expands the 
IL = o according to one of the two methods described in the 
previous sections. This is a good check of the procedure we 
adopted and is necessary when evaluating the first-crossing 
rate. 

The calculation of the first-crossing rate proceeds by 
integrating the probability density over S n and then taking 
the derivative with respect to S„. This is fortunate because 
we can directly compute 

2J / ddi . . . dS n didjdkWmb 

= ^w/l"^ n -° (90) 

We choose two different expansions for IT. The expansion in 
derivatives of Sect. 3.2 gives 

fB„ 



if 



dB3(S n ) 

2 



t(<>) 



r(«) 



2ttS 5 / 2 



d5n(n^ + ns + n^ +n^ + -- 



-ftr [-S + B 2 + SB' (B + 2SB')] 



-Erfc 



B 



B" 



(91) 



while the expansion using the approximation of Lam & 
Sheth (2009) (and discussed in Appendix A) gives 

fB„ 



If 



8B3(S n 



d5 n (u^ + ntt T) +ntr 



TVS* 



B (S n ) B„ 



On 



b>-, 



V(S n 



+ • 

V 2 (S; 

O n 



-B 2 (S„)/(2S„) 



where 



V{S) = J2 



{-S) p d p B(S) 



P ] - 



OSp 



(92) 



(93) 



Notice that the sum runs only up to p = 5 to provide a good 
fit to the data, as mentioned earlier in Sect. 3.1. If we now 
normalize the bispectrum as 



S 3 (S) = ^{8 3 (S)}, 



(94) 



we finally obtain the leading NG contribution to the first- 
crossing rate with a generic moving barrier. Using (91) we 



obtain 



■^ng(>S') 

Sz 

12 v / 2^S' 5 /2 



.-(o) 



.-(<•) 



T w + T w + T w + T 



(b) 



r(c) 



[-2 (S 2 + 2SB 2 - B 4 + SBB' (-75 + B 2 



-8S 3 B' 2 + 4S 3 BB' 3 



3V2nS 5 / 2 



S 



S 3 B"(B + 22SB')]e~ B /{2S) 
[B 2 + SBB' + S (-1 + 2SB' 2 )] e " s2/(2S) 

B 



-| ((2<S 3 + 55^) B" + SS 3 B'") Erfc 
while using (92) we obtain 

T- ^C^ B + V -B 2 /(2S) 



(95) 



+ - 



S, 



2ttS 3 / 2 
_ \B 4 -B 3 (P + 2SB') + 2B 2 (-S + P 2 

+SVB') + SB{V + 6SB' - 4P 2 B' - 2SV') 
-S(S + 2V{V + SB' - AST'))] e - B2/(2S) 



+ - 



3V2TS 5 / 2 



B 2 -BT- S 



-2P 2 ]e- B /< 2S » 



(96) 



where the prime denotes differentation with respect to S. 

Both formualae (95)-(96) can be further improved us- 
ing a saddle-point technique in order to resum the largest 
contributions from NG, as in D'Amico et al. (2010). Limit- 
ing this procedure to the leading terms of (96) and treating 
V(S) and the derivatives of -8(5*) as small parameters, we 
find for instance 



+ - 



■^ng(S') 

_P__ 

/2^S 3 / 2 6 
S 3 



Be 2s 1 

= — = e« 

V2^S 3 / 2 

-S 2 /(2S) 



S 3 B 



1SS 3 \ 
6 B J 



+ — —^ \-B 3 (P + 2SB') + 2B 2 (V 2 

+SVB') + SB(V + 6SB' - 4P 2 B' - 2SV') 
-2SV(P + SB' - ASV')] e - B2/(2S) 

B 2 /(2S) 



+ 



S 2 S 3 



\B 2 - BV - S + 2V 2 \ e 



(97) 

3V2nS 5 / 2 L J 

Notice that, in the limit of constant barrier, our formulae 
are slightly different from those of D'Amico et al. (2010); 
we believe that the origin of this difference is due to the fact 
that they assumed a very specific form for the cumulants 
(SiSjSk) oc (SiSjSk) 1 ■ With this assumption, one can find 
relations between the various derivatives of the cumulants, 
which otherwise are independent. 

In the limit of constant barrier B(S) — ^fa5 c one re- 
covers the spherical collapse result of MR3 (neglecting the 
terms proportional to £3) 



^{S) 



_Va6c_ e -ail/(2S) 

/2^S 3 /' 2 

SS 3 f(V^s c ) 4 



6y / a<5 c V S 2 



1 + 

JV^S C ) 2 
S 



(98) 



In Figure 5 we show the first-crossing rates (95) and (96), ap- 
plied to the case of the ellipsoidal barrier (2). The two curves 
differ by 0(10)% at most in the small halo mass regime. In 
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0.50 




5.(1 



Figure 5. The first-crossing rate deduced from Eqs. (95) (dashed 
blue line) and (96) (solid red line), for the case of ellipsoidal bar- 
rier (2). We used S3 given by Eq. (101) with local /nl = 100. 



1.08 



1.06 



1.02 




5.0 



Figure 6. The ratio between the non-Gaussian first-crossing rate 
/ng deduced from Eqs. (96) and the Gaussian one /q. We used 
cS3 given by Eq. (101) with local /nl = 100. 



Figure 6 we plot the ratio between the non-Gaussian first- 
crossing rate deduced from Eqs. (95) and the Gaussian one. 
In Figure 7 we show the ratios between the first-crossing rate 
given in (97) and the first-crossing rates (79) built up from 
two different commonly used form factors 7?.ng, the one of 
Matarrese et al. (2000): 



7?-NG 



cxp 




Uy/a6 e )S 3 d]ny/S 



(99) 



and the one of Lo Verde et al. (2008) : 

1 s r f(v^s c ) 4 



7?-NG 



1 + 



6 y/a8 c 



(yfa8 c 



\ S 2 



- 1 



+ 



dS 3 ({^8 c f 



dln-JS V S 



- 1 



(100) 



In the plots we used the conversion from the variable S to 
the variable M given in eq. (A2) of Neistein & Dekel (2008), 
while for the scale-dependence of 53 we used the following 
simple fitting formula 



S 3 (S) = 



2.4 x ID- 
S' - 45 



-,/n 



(101) 



which agrees well with LoVerde et al. (2008). 

As we can see, the first-crossing rate in the case of an 
ellipsoidal collapse and when NG is present is not generically 
given by the Gaussian first-crossing rate for the ellipsoidal 
model multiplied by the form factor obtained from the PS 
approach and can differ significantly from it by O(10 — 50)% 
or more at high redshift and large halo masses. 



5 CONCLUSIONS 



Excursion set theory provides an elegant analytical tech- 
nique to describe the distribution of dark matter in our uni- 
verse. When supplemented with various improvement con- 
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Figure 7. Ratio of the J-ng(S') in (97) to the first-crossing rate 
given by the J 7 st('5) ' n (51) times a form factor 7^ng> as a func- 
tion of the halo mass M for /nl = 100. The form factors are those 
in Eq. (99) (red lines) and Eq. (100) (blue lines). We considered 
redshifts z = 1 (solid lines) and z = 2 (dashed lines). 



cerning the physical modelisation of halo formation (such 
as the ellipsoidal barrier of Sheth & Tormen (1999) to take 
into account the triaxiality of halo collapse and the diffus- 
ing barrier of MR2 to take into account the stochasticity 
inherent to the process), as well as with improvements on 
some technical aspects (such as the inclusion of the non- 
Markovian dynamics introduced by the filter function), it 
provides a quantitative agreement with N-body simulations 
at the level of about 10% in most of the interesting mass 
range. While even more accurate results might be needed 
for precision cosmology, it is still remarkable that such a 
relatively simple theory catches quantitatively a significant 
part of the physics of such a complicated dynamical process 
as the formation of dark matter halos. The same is true if 
the excursion set method is applied to describe the abun- 
dances of cosmic sheets and filaments. In this paper we have 
extended the path integral approach proposed in MR1 for 
the spherical collapse case to the case of generic moving bar- 
riers using a top hat window function in wavenumber space 
. We have shown that, using a well controlled and system- 
atic expansion, we can reproduce the ST halo mass function 
very well, therefore putting it on firmer grounds. We have 
also performed the computation of the first-crossing rate for 
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the ellipsoidal barrier in the presence of non-Gaussian initial 
conditions. Our result is given in eq. (97): it is fully consis- 
tent in the sense that it does not require the introduction 
of any form factor artificially obtained from the PS formal- 
ism based on the spherical collapse and in fact it provides 
a halo mass function which quantitatively differs from the 
one obtained from the form factor procedure. 
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APPENDIX A: REPRODUCING THE 
FIRST-CROSSING RATE OF SHETH & 
TORMEN 

We first compute n (1) . Using eqs. (24), (29) and (30), the 
expression of n^ (8 n ; S„) in eq. (48) can be rewritten as 



r(!) 



U\t! (5 n ;S n ) = 



B n (B n - *„) v (-l) p 



7T 



E 



n\ 



-B 



(p) 



(Al) 



P =i 

S„ (q c.\p-(3/2) 

>0 



S 



3/2 



xe -Bi/(2S i ) e -(B n -S n ) 2 /l2(S n -S i )] _ 

Instead of computing directly this integral, we now re- 
call that to compute the first-crossing rate (39) we need 
to compute the first derivative of n e (<5„; S„) evaluated at 
8 n = B(S n ). Since the integral in eq. (Al) is finite in the 
limit 8 n — > B(S n ), taking the approximation (S n — SiY^ 1 ~ 
[Sn) v ~ x does not alter the convergence properties of the 
integral, but simplifies significantly its computation. This 
is equivalent to the approximation made by Lam & Sheth 
(2009), see in particular the discussion below their eq. (20). 
Exploiting the fact that 



dS, 



sf /2 (s„-s0 1/2 



xe -B J (S„)/(2S s ) e -(S(S„)-i„)V(2(S„-5 i )) 



^/2^ : 1 
BJSZ)£T^ 



(25(5; 



n I u n ! 



Zo n 



(A2) 



The former does not contain a finite part in the continuum 
limit and its divergence cancels against the divergent part 
of the latter sum (see appendix B of MR1). Thus, we are 
reduced to compute the finite part of the sum over i < j. 



C 1 )/ 



Proceeding as before for the calculation of IT (S„; S„), and 
taking again (S„ — Si) 11-1 ~ S^ 1 we obtain 



r(2,ST) 
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(A4) 



Let us indicate the integral by A(S n , S n )- It is convenient to 
perform the inner integral by deriving with respect to 8 n 
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(A5) 



where we used eq. (B.26) of MR1 in the second line and 
eqs. (A. 5) of MR1 and (A2) in the third line. Integrating 
over 5 n we find 



A(Sn,S n ) = 



2n (B(S„) — Sn) -(2B(S„)-.5„) 2 /(2S„) 



Si' 2 B(S n ) 
which can then be inserted into eq. (A4) to give 
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contributions II 



(n,ST) 

<! = 



we find that 11^ ' (S„; S„) (where the superscript reminds A similar procedure can be used to show that higher order 

us that we have approximated the integral) is given by 

II (1 _: ST> (8 n ;S n ) = 2 (B{Sn) - S n ) ^(2B(S„)-<S„) 2 /(2S„) 



2-kSI' 2 

E (- S ™) P R (p) 
p=l 



(A3) 



(n > 2) vanish as (B(S„) — 8 n ) n when 
8 n approches the barrier value B(S n ). 

The calculation of the first-crossing rate is then straight- 
forward, through eq. (39). The zero-th order contribution 
from n^Jg is given by eq. (68), while the first-order contri- 

(1 ST) 

bution from II^_! reads 



Next, we compute IT {5 n \ S n )- The sum over i,j in eq. (49) 
can be split into a sum over i = j and a sum over i < j. 
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Q£ (B n — S n ) 2 



Higher-order contributions to the first-crossing rate vanish 
This is already clear from the contribution arising from the 
second-order IT 
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r 



(2, ST) 



(Sn) 
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which vanishes for 8 n = B„. The total first-crossing rate for 
a moving barrier, in the approximation discussed above, is 
therefore given by 



Jst(5) = 
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APPENDIX B: COMPUTATION OF IT^,, 11^,, 
n (c) 

In this appendix we compute the contribution to n t= o in 
the derivative expansion discussed in Section 3.2. The first, 
using the techniques discussed in MRI, is simply computed, 

n<?2 (5 n ;S n ) = --^B n (B n -S n ) 
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2\ 1 ' 2 dB n (B n -S n ) 
' ■ exp 



^' dS n Si' 2 
The second term is 

I d 2 B 



B n (B n — S n ) 

Abi Zybn bij 

(2B n -S n ) 2 



Zo n 



(Bl) 



ni%(s n ;s n ) 





I d 2 B. 
2n~dSl 



-.3/2 



2 \B n — 5 n ) 



2rr dsi Bn(Bn ~ Sn) 



B n 
'2Si 



, (Sn bi) 

dbi ^ exp 



(Bn - Sn) 2 

Z\u n Ji) 



(B2) 



2nS 1 J 2 e^ 2B -- s - ) /(2S "» 



Tr_B n Erfc 



/ 2B n — Sn 



' 2b n 



where the integral has been computed using eq. (109) of 
MRI. The last term is the most complicated. Using the a- 
regularization and the finite part prescription developed in 
Appenix. B of MRI, we find as usual that the terms in the 
sum with i — j have a vanishing finite part, while the contri- 
bution from the terms with i < j (plus an equal contribution 
from i > j) can be written as 

2 



r(c) 



n^ (*„;S„) = 



B n (B n — 8 n ) / 



dBn 



7rV27T 

xTP dS t / dS 

Jo J Si 

Bn 



x exp • 
_ B n (B„ 



V dS„ 

\bn bi) 
^(^-SiF/HSn-S,)*/* 

ae {B n — 5n^ 



2 rS, 



x exp 



TT^/2ll 

BZ, 



S_n)_ (dBn 

\d,S n 



dS, 



\b n bi) 







s 



3/2 



(B3) 



-^ I ^ J*. ^ {Sj - Si)W(Sn - 5,)V2 



x exp 



Z\Oj Oi) ZyO n &j) 



(B4) 



where TV denotes the finite-part prescription developed in 
App. B of MRI. The integral over dSj is performed using 
MRI, eq. (108), and is equal to 



2tt 



Se (Sn - Si) 1 ' 2 



exp 



(B n - S n + y/ae) 

2[o n bi) 



(B5) 



Expanding the exponential we therefore get a singularity 
l/-\/i (which is canceled by a similar singularity in the term 
of the sum with i = j, see MRI), and a finite part, given by 



'2n 



(B n ~ Sn) 

(S„ - Si)*/ 2 



--(Bn-Snf/mSn-Si)] 



(B6) 



The remaining integral over dSi is performed again using 
MRI, eq. (108), so finally 



r(c) 



U.^ (S n ;S n ) = 



2\ 1 / 2 _ „ ,2 (dB r , 

\dS n 



(B n — 8 n 



C.l/2 

O n 



x exp 



{2B n — 5<n 
Zo n 



(B7) 



APPENDIX C: COMPUTATION OF n[% 

In this appendix we fill the missing step in the computation 
of n);J . The issue is the computation of the integral 



Xp(a,b,S n ) 



x exp 



dbi b^ \bn bi 

b 2 



a 
'2St 



Zyj n Oi) 



(ci) 



where a = B n > and b = (B n — S n ) > 0. Changing the 
integration variable to z — (S n /Si) — 1 we get 



lp(a,b,S„) = S n exp < — 



a 2 + b 2 

Abn 



dz 



1 1 

2 3/2 z l/2 



1 + Z 



x exp 



a 

2b r , 



Zjn 



(C2) 



For p — 0, 1 the integral can be performed exactly (see 
cq. 9.471.12 of Gradstein & Ryzhik (1980)) and we get 3 



lo(a,b,S n ) 
Xi(a,b,S n ) 



( 2?r ) a + b -(a+6) 2 /(2S„) 



s; 

(2rr 



?,/2 



,V2 



ab 
1 _ 



(a+6) 2 /(2S„) 



s, 



1/2 



(C3) 
(C4) 



For p ^ 2 we have not been able to compute the integral 
exactly. However, as discussed in the text, for computing the 
first-crossing rate it is sufficient to evaluate it at b — 0. The 
resulting integral can be computed (e.g. using Mathematica) 



3 These integrals were already computed exactly in a different 
way in MRI. We thank Ruth Durrer for suggesting this more 
direct derivation. 
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in terms of the confluent hypergeometric function U(a, b, z), 



l p (a,0,S n ) = S p n - 2 ^e- a2 ^> 



xTIp- -)U[p-l 



1 a*_ 



(C5) 



Observe that U(0,b,z) = 1 and T(l/2) = 0F, so cq. (C5) 
also reproduces correctly T p (a, 0, S n ) when p — 1. It is also 
useful the limit 



J r PlimX p (a,0,S„) = -7rc p ST , (C6) 

a— >0 



X p (0,0,S n ) 

where the coefficients c p are given by 

2 r( P -i) 



APPENDIX D: COMPUTATION OF THE 
GENERAL TERM U ( ^ IN THE LIMIT 

(B N -S N )^0 

The general term ff^g is given by 

! ~ B (Pl)... B CPm) 



n 



(m) 



E 



n— 1 



,...,i m =l 

d8 1 ---d8 n - 1 d ll ---di m W sm . 



The last integral is equal to 

f Bn dSi---dS n -idi 1 ---Qi m W m = 

J -co ± m 

n gm (So , B n , s n )n gm (b„ , B n , s l2 - s n ) ■ • ■ 

• • • n gm (B„, B ni 5 im - S lm _ 1 )n gm ( J B„, «5„, S„ - 5,„ 
Using eqs. (29)-(31) for n gm , cq. (Df ) becomes 

CO 
(m) J- J^n\J^n On 



ll 



m! 2^7^ 



E 



\pi+...+j>„ 



Pl,-,Pm=l 



w 



r(pi) ... R 



+C(B„ - 8 n f , 



<TVm) (T3 o 



where 



r (m) 



Jo 

l-S n 



S n 



(C7) 



(Dl) 



(D2) 






X 



■ />',„ ^" T J ,.,,. . (D3) 



(&, 



-£*, 



")3/2 



We have only considered the finite parts from the sum 
with ii < 12 < • • • < im, because the divergent parts all can- 
cel. A priori, we cannot exclude that there may be other 
finite contributions to the sum coming from terms with 



ii < • • • < ik — ik+i < ' ' ' im- However, we expect the con- 
tribution we compute here as representative of the correct 
result. 

The integral Jp\}... Vm (B„, S„) satisfies the recursion re- 
lation 



r (m) 



7 ym) (R <? 
'Jpi ,...,p m V n ' tJn 



Let us set 

r (m) 



s n 



dS, 



; 3/2 



o S t 

'"■ ' m w - <-. i (D4) 



^ ^P2>---jPm \ ' n *) ' 



Jp™'.., Pm {o,y) = (-7r) m c P1 ,..., Pm V 



pi+.-.+Pr, 



(D5) 



where the coefficients c are now to be determined. We in- 
sert the ansatz above into the recursion relation (D4) for 
Jp™~Jp m {Brn S„) and obtain 

>^Plv,Pm+l (Bn, bn) = ( — 1") c P2,---,P m +l 

x / S " dSi(s n - s^ + - + ^-^s; i/2 e~ B ~'^Xm) 

The previous integral is solved with the substitution z — 
(Sn/Si) — 1 and it evaluates to 



l+l TT. + 3 

Pk o— 



v^™+i _ m _ 9 Y^ m +i 

nZjfcl Pfc 2 roc , z £*ik=l 



-(1+^) 



— &n 






2 / 



; 2S : r Cft 



therefore eq. (D6) becomes 



, ,, -,_,,,. , fc=i PA; 2 

^PlvjPm+l V-"^! °"7 — V 'V L -P2,---,Pm + l Jn 



^ A / ^^ e 



^(Er> fc 



2 , 



XC/ 



/Y^tti + I 



Pk 



n+i i _b; 

2 ' 2 ' 2S r , 



(D8) 
0, and 



We can evaulate eq. (D8) in the limit B 2 l /(2S n 

retain the finite part only (as the divergent terms all cancel 

in the end): 



J^l]L +1 (o,v) = -2v^(-^ m Cp2 ,..., Pm+1 
r (E fc=1 ~-t-i) 



(D9) 



On the other hand, the left-hand side of the previous 
relation can be expressed by (D5) and we then arrive at 
a recursion relation for the coefficients c (after relabelling 
m — »• m — 1 for convenience): 



2 r(Er=i^- 

lid for r 
37) 

2 r( P -i) 



m+l ' 

2 , 



L -P2,---,Pm I 



(D10) 



which is valid for m ^ 2, while for m = 1 we have already 
found in (C7) 



(Dll) 



Vi r (p - 1) 



Equations (DlO)-(Dll) define recursively the coefficients c 
and it is possible to find them easily up to any desired order. 
As the c appear in the generic integral (D8), which in turn 
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appears in (D2), it is then possible to write down the result 
for the generic term FT" 1 ' : 



n (m) (B n - S„)e 2S ' 

11 e=0 — ~ , m _, 3-m 

ml 2 2 'i ! 

d(pi) r>(Pm) 



oo 

y" (-i)Er=i pfc+m+i 



oi,...,p m = l 






x r 



E 



p* 



t/ (X>- 



m + 1 1 B;J 

.Z .Z ZnJn 



+ 0(B„ - 6„ 



(D12) 
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